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1. Introduction 

To study the nature of matter in QCD at high temperature T or large quark chemical potential /I 
is one of the most challenging issues in particle physics. Several novel phases have been proposed, 
such as quark-gluon plasma and color superconductivity. Precise determination of the QCD phase 
diagram on the (/I , T) plane will provide valuable information for experimental search. 

In the continuum, the grand canonical partition function of QCD at finite T and jl is 



where N is particle number operator N = J d 3 x (x)\}f(x). 

Monte Carlo (MC) simulation of lattice gauge theory (LGT) is the most popular nonpertur- 
bative method based on the first principles. This approach has been successfully applied to QCD 
at finite T with zero jJ.. However, LGT experiences serious problems, like species doubling with 
naive fermions and complex action at finite jU. 

The Hamiltonian formulation of LGT at finite pL does not have the complex action problem^ 
||[ ||, |J. The complex action problem in Lagrangian formulation forbids numerical simulation at 
real fx. The recent years have seen enormous efforts [||, [6|, [7], ||] on solving the complex action 
problem. 

There have been several popular approaches to solving the species doubling problem of naive 
fermions. The staggered (KS) fermion approach preserves the remnant of chiral symmetry, but it 
breaks the flavor symmetry and doesn't completely solve the species doubling problem. The Wilson 
fermion approach avoids the doublers and preserve the flavor symmetry, but it explicitly breaks the 
chiral symmetry; In order to define the chiral limit, one has to do nonperturbative fine-tuning of the 
bare fermion mass. 

The overlap fermion approach^, [RJ is claimed to have the properties that chiral symmetry is 
preserved and species doubling problem may be solved. However, the Dirac operator is nonlocal, 
and the computational costs for simulating dynamical overlap fermions are typically two orders of 
magnitude heavier than for the Wilson or KS formulations. It is also very tough to introduce the 
chemical potential into the action. Before the breakthrough of numerical algorithms for applying 
overlap fermions to QCD thermodynamics, it is very useful to do an analytical study. 

In this paper, we summarize our study on above issues using Hamiltonian LGT with Wilson 
fermions, and Lagrangian MC simulations with four flavors of Wilson fermions. We also present 
some new results from strong coupling analysis of Lagrangian LGT with overlap fermions. 

2. Hamiltonian lattice QCD with Wilson fermions 

We begin with the QCD Hamiltonian with Wilson fermions H = H g + Hf at jU = on 1 
dimensional continuum time and d = 3 dimensional spatial discretized lattice, where 




(1.1) 
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(2.1) 



According to Eq. (1.1), the role of the Hamiltonian at finite jU is played by = H — fj,N. 
Denoting Nf, N c , V the number of flavors, colors, and spatial lattice sites, and |£1) the vacuum state 
when (Q.\H^\Q.) is minimized. For free massless Wilson fermions at T = 0, we obtained[|l]] the 
energy (£l\H = 1N c NfY, p \p\ [© (M ~~ \p\) ~ 1]> an d the subtracted energy density in the infinite 
volume limit V — > °° and continuum limit a — > 



M=0 



N c NfV 



(27T) 



871 

(2^ 



3, 

" ^=4^' 



(2.2) 

which agrees with the continuum theory. Therefore in the lattice Hamiltonian formulation, the 
chemical potential could be introduced in a natural way as in the continuum. 

For infinitely strongly interacting Wilson fermions, integrating out the gauge fields leads to 
four fermion interactions^]. Extreme conditions (large T or jJ.) induce chiral phase transitions. For 
Nf/N c < 1 with N c = 3, we obtained an equation for the critical line where the chiral condensate 
and the dynamical mass of quark vanish continuously^] 



(1+r 2 ) 



2T' 1 + 
c +7 ^ ln 



27 c 

l+3r 2 



l + 3r 2 



1 



27^ 

l+3r 2 



(2.3) 



Here we have rescaled the chemical potential and temperature as \i' = fl/(3K/a) and T' = T / (3K /a), 
with K being the effective coupling of four fermion interactions. Below some r 3 ', there is a first 
order chiral phase transition line[^] 



^=l+2r 2 



(2.4) 



where the chiral condensate and the dynamical mass of quark vanish discontinuously. The point 



when two lines described by Eq. (2.3) and Eq. (2.4) join at lower T 1 is the tricritical point, as 
shown by the filled circle in Fig. [I]. 



3. Lagrangian lattice QCD with Wilson fermions 

The lattice action at }X = proposed by Wilson[[TI|] is S = S g + Sf, where 

S g = -^Tr(U p + U;-2), 
p 

s f = LvMMjco-v'GO, 



.Y.V 



M, 



.Y.V 



(r-Yj)Uj( X )8 XA ,_j+(r + Y j)ul(x-j)8- 

7=1 L 



x,y+j 



(3.1) 



with /3 = 6/g 2 and K = 1 / (2ma + 8r) . The lattice Hamiltonian H in Eq. (2.1) could also be derived 
from the Wilson action by Legendge transformation. 
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Figure 1: Phase diagram from Hamiltonian lattice QCD with massless Wilson fermions at strong coupling. 
The solid and dotted lines stand respectively for the first and second order transitions. The circle is the 
tricritical point. 



However, naive introduction of the chemical potential would lead to divergent e su b in the con- 



tinuum limit. In Ref. [|12j, the chemical potential is introduced by replacing the link variables in the 
temporal direction in fermion action in Eq. ( |3.l| ) with L^*) — ► e a ^U^(x) and ul(x) — ► e^ a ^lll(x). 
The fermionic action is reduced to the continuum one when a — > 0. 

Nevertheless, the effective fermionic action in the partition function becomes complex, and 
forbids MC simulation with importance sampling. Several revised methods, e.g., improved reweight- 
ing [j|] and imaginary chemical potential methods^], 0], were proposed to simulate QCD with KS 
fermions at finite /x. 

Lattice QCD at imaginary chemical potential does not suffer the complex action problem. 
In Ref. we applied this method to the MC study of the phase structure of Nf = 4 Wilson 
fermions. We measured the expectation of the Polyakov loop, chiral condensate and their suscep- 
tibilities for various (Hi,T) at some K. From the position of the peak in the susceptibilities, we 
determine the transition point. Replacing fj,j by — ifi, we directly continue the transition line on the 
{Hi,T) plane to the real (ju, T) plane. 

Figure ^| is the expected phase diagram of lattice QCD with Wilson fermions in the (ju , T, k) 
parameter space. There is a surface K = K c hi ra i where the pion becomes massless. Above this 
surface, there is no phase transition, as confirmed by our numerical simulations for K = 0.25. 
Interesting physics is below this surface. Of course, the order of transition depends on the value of 
K. For r = 1, we find JQ G (0.001,0.15), fc 2 G (0.15,0.165) and KMnd G (0.17,0.25). 



4. Lagrangian lattice QCD with overlap fermions 

The overlap fermionic action Sf has the form|jl(J 

S f = m£v/(x)vA(jc) + ^^(jc)D(x,y)vA(y), 



151/4 



Phase structure of lattice QCD with Wilson and Neuberger quarks . 



Xiang-Qian Luo 




Figure 2: Phase diagram of lattice QCD with four flavors of Wilson quarks in the (jU, T, K) parameter space, 
from MC simulations. For KG [0,K\] and K £ [Kz, Kchuai], the phase transition is of first order, and while for 
K E (fCi , K2), the transition is a crossover. 

D = 1+X^L=, (4.1) 
where a is set to be 1 for convenience. The operator D is nonlocal and it is extremely difficult to 



do analytical calculations. In Ref. []15[], we used the Taylor expansion trick[]13|, |14| ] to derive an 
overlap action at finite jit 

m\ C J, 



Sf = (l + ^)^(j^Z[q(x)rjUj(x)q(x+d 

+ Y}e^q{x)y A U 4 {x)q{x + A)-e-^q{x + A)YAUl{x)q{x)]\ +m^q{x)q(x), (4.2) 

X J X 

where A = 4r — Mo, and C = t/2, with t an expansion parameter. The fermion fields q and q are 



related to y and i// by [14] 



q = xjr, g = M-IdJ y. (4.3) 

The chiral order parameter is then given by 

(qq) = (y(l-\D) ¥ ). (4.4) 

In Ref. jl5|], we studied the phase structure of LGT with overlap fermions on the (jJ,, T) plane 
at the strong coupling. The phase diagram is shown in Fig. ||. We find that the phase structure is 
very similar to the Wilson fermion case. 

5. Discussions 

LGT with imaginary chemical potential could be used for simulating QCD at small pL. The 
overlap fermion approach is a promising one for investigating chiral properties of the phase transi- 
tion, and the Hamiltonian lattice formulation is a more natural way to introduce chemical potential; 
However, a lot of efforts have to be made before realistic MC simulations could be carried out. 
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Figure 3: Phase diagram of lattice QCD with massless overlap fermions at strong coupling on the 
plane. The dotted and solid lines stand respectively for the second and first order transitions. The circle is 
the tricritical point. 
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